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УДК 631.354-530 


Принцип минимакса в теории и практике оптимального проектирования 
приводов современных зернокомбайнов 


В. В. Радин 
(Донской государственный технический университет) 


Поставлена и решена задача оптимального выбора конструктивных параметров элементов приводов зЕРно- 
уборочных комбайнов, когда эти элементы реализуются с ПОМОЩЬЮ клиноремённых и гидравлических пере- 
дач. Для этого использован принцип минимакса, реализованный с помощью метода оптимального проекти- 
рования параметров динамических систем в пространстве состояний под действием переменных нагрузок. 
Для формирования пространства состояний использована ранее полученная автором теоретическая модель 
динамической системы привода комбайна с нелинейными неголономными связями первого порядка. 
Ключевые слова: динамика привода зернокомбайна, нелинейная неголономная связь первого порядка, 
метод наискорейшег О спуска, оптимальное проектирование. 


Введение. СА|5-технологии, стремительно развивающиеся в мировом комбайностроении, требу- 
ют так выбирать параметры отдельных элементов приводов технологических машин, чтобы обес- 
печивался минимум максимальных динамических нагрузок на рабочих органах. Предлагаемые 
ниже метод и алгоритм оптимизации параметров приводов зерноуборочных комбайнов на основе 
принципа минимакса являются актуальными и могут быть полезными для соответствующих конст- 
рукторских организаций. 

Постановка задачи. Рассматривается четырёхмассовая динамическая модель привода зерно- 


уборочного комбайна (рис.) с клиноремёнными и гидравлическими передачами. Приведённые 
функции передач В,, й =1,3, представлены в данной модели как нелинейные неголономные свя- 


зи первого порядка. Уравнения основного движения четырёхмассовой системы построены по ме- 
тодике, представленной в работе [1] и применимой к сколь угодно сложному приводу зерноубо- 
рочного комбайна. 

Пространство состояний динамической системы, адекватное реальной нагруженности про- 


ектируемого привода, описывается реакциями В; (2) неголономных связей. Пусть В, (2) — век- 
1[0;№] 


тор-функции, заданные и непрерывно дифференцируемые на некотором открытом множестве 
СЕ, 

где Е, — евклидово пространство. 

Предположим также, что нам задано выпуклое замкнутое ограниченное множество 

со’, 

представляющее собой максимальные динамические нагрузки на элементах кинематических 

трансляторов проектируемого или модернизируемого комбайна. Уравнения основного движения 

системы выделенных роторных масс этого комбайна, изображённые на рисунке, запишутся в нор- 


мальной форме уравнений Коши следующим образом: 
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Ф-РЕГ, (1) 


где 2(Ё) = [Ф, (К Ф.(8 Ф.И ©. (#)] — вектор переменных состояния динамической системы 


(см. рис.); Т — символ транспонирования вектора состояния; Т» — время наблюдения; 6 — 





вектор переменных проектирования, например Ь = [ А, 5, Е, Г] ‚ ВЕН+ШП; ГЕЪП; 


п = 4; р = 1; А» 5, Е, — параметры клиноремённых и гидравлических передач приводов про- 
ектируемого комбайна, определяющие невосполнимое отставание ведомых элементов от ве- 


дущих в А-м контуре; 1;— приведённые к соответствующим валам моменты инерции выделен- 


ных в динамической схеме масс (см. рис.) относительно их осей вращения; (2,2 (Е),Б) И 


вектор правых частей дифференциальных уравнений (1); 2’ — вектор начальных условий. 


К. 1 


|, С. 
мт ФЕ Же © 


м ГЕ МУЗЕМИ| М5 


ооо 
. Ф : Ф 


Динамическая модель привода комбайна: ф,, ф,, ф,, ф. — угловые скорости соответственно двигателя, ведущих колёс 
































комбайна, главного контрпривода, молотильного барабана; Г, — приведённые моменты инерции соответствующих эле- 
ментов приводов; М, — приведённые обобщённые силы на соответствующих валах; В, — приведённые функции передач 


как связей 


Полагаем теперь, что самым энергоёмким рабочим органом зерноуборочного комбайна 


является молотильный барабан, и обозначаем через ф, его оптимальную угловую скорость, 


обеспечивающую наилучшее качество выполнения техпроцесса. 
Задача минимакса будет состоять в том, чтобы вектор Ь из пространства переменных про- 


ектирования был выбран таким образом, чтобы обеспечить минимум максимальных вариаций: 
глах |ф, (Е,5) -Ф‚|= п, = [0,Т, |, (2) 


где ф, (Е.Б) — текущая угловая скорость молотильного барабана. 
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Если принять ® (Ё,6,2) = Ф.-Ф.| и $, =тах® (ЁЕЪ,2) в момент времени Е е [0,Т, |, то 
задача (2) конкретизируется: 
плитах ® (#,6",2) = миф, Ё=[О,Т, |, (3) 
Вычисление тах ®, (Е,Ъ,2) заменяется точной верхней границей ф, = В,,, которая опреде- 
ляется соотношением 
$. (-Ф 56, Ее [0,Т, |. (4) 
Тогда задача (2) сформулируется следующим образом. 
Требуется определить 





тиф, (5) 

при ограничениях 
| (Е,6) -Ф.|< В, Ее [0,Т, | (6) 
< <Ь’, 1=1113, (7) 


где 6’, 5” — наименьшее и наибольшее значение переменных вектора проектирования. 


Для данного проекта в соответствии со схемой (см. рис.) и формулой (1) максимальное 
число переменных проектирования равно 13. 

Для решения поставленной задачи (5) с ограничениями (6) и (7) используем методику оп- 
тимального проектирования параметров динамических систем в пространстве состояний под дей- 
ствием нестационарных нагрузок [2]. 


Заменяем ограничение (6) интегральным ограничением вида: 
И 
ф= | (© (6.2.6) -6.)@ =0, (8) 
0 


где подынтегральная функция определяется соотношением 


#(Е,2,6)-В., Е (Е,2,6)-Б. >0 
ель -в,)- | о 


0, если Ф, (Е,2,6) -Ь,, < 0 и 


Если интегральное ограничение (8) в ходе реализации метода наискорейшего спуска не 


выполняется, вводится в рассмотрение сопряжённая система дифференциальных уравнений: 


т 


Е | [9] 
А+. А- эп ((^, — 6. +=) 


, (10) 
А (Т,) =0 
где ^(Ё) — вектор сопряжённых переменных, равный 
(Е) = [Л, (2) ^, (Е) А: (2) ^. (#7, (11) 
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№ { 
= — матрица Якоби для вектора правых частей системы дифференциальных уравнений 
7 





(1); 7 =[0 о 0 $9п(Ф, (Е) - Ф, (| ‚ так как ^ =. (Е, 2,6) -Фо|; = — =-активные огра- 


ничения, вводимые для получения устойчивых коэффициентов чувствительности, которые 
вычисляются, когда в вектор проектирования включена и точная верхняя граница В. 
— т 
ь 
Ь=| 5’ В, ]|, 


оф 
0Ь,, 





а таблица коэффициентов чувствительности будет включать производную 


т 


И 
= -| Яап - В+) а |. (12) 
0 


Конкретная система дифференциальных уравнений основного движения (см. рис.) пред- 
ставляется следующим образом: 
ф, =2 (6,2,Ё) = [Ми (6,2,6) + У, (,2,6) 9, (6,2,Ё)+ У; (6,2, Е) 9, (6,2,Ё) + 
+Ун (6,2, Е) 9 (6,2, ) | /|1— М,» (6,2,Ё) 9.1 (Ь,2,Ё)- М3 (6,2,Ё) 9,» (Ь,2,Ё)- 
Уз (6,2,) 9, (6,2, Ё) - 9., (6,2, | 13) 
Ф, = 9. (6, 2,6) (Ь,2,Ё)+ 9, (6,2,Ё) 
ф, =9,, (Ь,2,Ё)Е (Ъ,2,Ё) + 9,, (6, 2,Ё) 
ф, = 9.1 (62,6) Е (Ь,2,Ё) + 9. (Ь,2,Ё) 


В системе нелинейных дифференциальных уравнений (13) принято: 








С: Л.Е,Ф, 
ав. И: 
11 2 12 [о 
№ С Са —С: Са : _ 92Сы Г 93С» : 
92: = АЛ ' 22: АЛ ' 
= С Сзз —С:Сз . ‚5 92Сзз = 92Сзз : 
31 АЛ ' 32 АЛ ' 


ААС Сивас 
м = (м, + С,.б,› +С;6., +С..0. )/т, з 
а и: 
Моно ИТ 


9, =-Лз ГА, + Л (1-2=. )^, |ф,; 
Чз =-Л |(1-=, А [Е, + Л (1-2=, 2, + (1-=.) АР]; 


5. са ие: 5, = М, /С,.; 5. = Г. С„/ МЛ; 
б. =-1С.з /АА; 6.1 — С» / АА; б-› = [4С» /АЛ; 
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5. = (М.С; — М.С )/А^; 5.3 2 (М.С; — М.С )/АХ; 


С; — элементы матрицы С размером 3х4: 


ий (1 =” =) Е 724В,ф, 0 0 
С = —Лз (1 = =) 0 т 2) :%В-Ф.Ф: 2 Л ЛАВ (1 и 2=.)ф.Ф. 
Ла (1 Е =. ) 0 2).АВ,Ф.Ф. (1 `. =.) т 2) .ВФ!Ф. [^, (1 — Е) ы Л.А (1 — 24 ) (1 № 2=. ) 


к — Расчётные передаточные отношения между осями валов /го и К-го элементов модели (см. 
рис.); =, — функции невосполнимого отставания А-го элемента модели от /го. 
В настоящем исследовании принята линейная зависимость функций =, от нагрузки М, на 
[м валу: =, =АМ,. 
Крутящие моменты на валах привода комбайна, входящие в систему (13), определяются 
выражениями: 
М, =Е, + Е, + В,Ф,; 
М, = Е, + Е, + ВФ; 
М; =Е, + Ем, + В.Ф), 
где Г, — составляющие крутящих моментов от сил полезного сопротивления; Г, — составляю- 


щие крутящих моментов от сил вредного сопротивления; В, — коэффициенты, определяю- 


щие составляющие крутящих моментов на валах, пропорциональные квадратам угловых ско- 
ростей этих валов. 


В свою очередь, составляющие крутящих моментов РЕ, распределены в пространстве со- 
стояния рассматриваемой системы по фазе воздействия, которую можно оценить величиной чис- 
того (транспортного) запаздывания т, следующими соотношениями: 

Е =С+К,а (Ё-т,); 
4 
. = 2.Кы9 (Е № р 
Е, =К.9(Е-т.); 
4 (Е) — подача хлебной массы, С = соп&. 


Невыделенные элементы привода в данной модели приведены к валу главного контрпри- 
вода. Их число может быть любым. Однако адекватные пространства состояний привода комбай- 
на генерируются восьмимассовой обобщённой динамической моделью [1]. Поэтому невыделенные 


в схеме элементы привода принимают значения: 


пе а. 


Коэффициенты К’, / =2, 3, 4, определяются соотношением: 
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—_ УМ, „К 


К. 
94° +Ко4 


/ 


где у, — коэффициенты, характеризующие распределение внешних нагрузок по Ам компонентам 


схемы (см. рис.); 4” — расчётное значение подачи хлебной массы, определяемое техзадани- 


ем на проектирование; М, — крутящий момент двигателя в точке стыка регуляторной и кор- 
ректорной ветвей внешней скоростной характеристики; К» К, — моделирующие коэффициен- 
ты; 4 — начальная средняя величина подачи хлебной массы. 

Крутящий момент М, развиваемый двигателем комбайна, задаётся аппроксимацией внеш- 


ней скоростной характеристики: 
2 , 
ыы Е 
5». а, Фит 7 Ф, > Ф!: ти 
1=0 


2 
о Е 
М, = У аф, Фи < Ф, <Ф ; 
1=0 
а(Ф. -ф,), Ф: = < Ф, < ЧФ: тах 
0, Ф, > Фатх 
где ф,, — угловая скорость двигателя в точке стыка линейной и нелинейной частей внешней 


скоростной характеристики; а; — коэффициенты аппроксимирующего многочлена нелинейной 
части внешней скоростной характеристики, определяемые методом наименьших квадратов; 
59 а — коэффициент статизма регулятора угловой скорости двигателя, задаваемый при выбо- 
ре мощности двигателя в соответствии с техзаданием. 
Уравнения нелинейных неголономных связей кинематических трансляторов в проекте 
представлены в следующей форме: 

Ф, — Л> (1 = =, )Ф, 

Ф: — Лз (1 -=з)Ф, 

Ф. > Л. (1 - =. ) (1 - =. )Ф, 


ОЕ . 
Матрица Якоби 7 вычисляется как производная сложной функции со многими перемен- 
2 


ными по следующей схеме: 


М, 
9С 58 , ‚ 4 ‚ ЗА м ‚ 56 м ты (14) 
092 092 097 092 092 02 02 02 





д 
Далее вычисляем матрицы 56 . При вычислении этих матриц используется та же последо- 


вательность, что и в схеме (14). 


В проекте принято: 


м, бе 
=-.М. —@; 
ааа РУ 





дм, _ ОМ, _0М. _ а 
06 06 06 ' 0 
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0, мА, ам, 
ЕТ. 06 ``” 0Ь 





1=2,3, 4. 


Ь я [, Ь, Ь, Ь, Ь, Ь, Ь, | Ь. Во Ь > Ь ] ' 








81 0 ообобооббо0тоо0", 
6 -" - 
91 о обобооббоботоо", 
95 - 
913 о обооббобоббот0", 
РИ - 
1 о обооббобб0б0б001т. 
РИ : 
Выводы. 
1. Метод наискорейшего спуска обеспечивает устойчивую сходимость решения по- 


ставленной задачи при использовании диагональной весовой матрицы переменных проектирова- 


ния Г,, / = 10,13, следующего типа: 


И =[0 20202047. 


2. Применение изложенной выше методики к оптимизации конструктивных парамет- 
ров элементов привода зерноуборочного комбайна позволило установить, что минимальные ди- 
намические нагрузки в приводе достигаются при выборе следующих технических параметров: 

— мощность двигателя комбайна — 250 л. с. 

— приведённый к валу двигателя момент инерции привода комбайна с измельчителем со- 


ломы — Г, =15,7 кг.м. 


— коэффициент статизма внешней скоростной характеристики двигателя 35,3. 
— резонансные для основного движения частоты внешних возмущений со стороны подачи 
хлебной массы находятся в диапазоне 0,85—1,3 Гц. 
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МТАМТМАХ РЕТМСТРЕЕ 1№ ТНЕОКУ АМО РВАСТТСЕ ОЕ МОБЕВМ СКАТМ СОМВТМЕ АСТОАТОК 
ОРТТМАЕ ОЕ$ЗТСМ 


\. У. Вадт 
(Роп За ТесНтса! Упмег®Ку) 


Тре рго Мет оп {те деяюп иа!иез орйта! спосе оЁ пе дгат сотЫте аМие еетепё, ипеп езе еетепё5 аге 
/трететеа {7гоидв пе У-БеЁ апа пуагаш!/с {гапзт/$!опо, [5 5еЁ апа зо/уеа. Рог #5 ригрозе те титтах рип- 
сре, ‘трететеа ргоидйв пе те{тоа о {те орита! ае$!дп оЁ {те зЕйе-зрасе дупат!с зузет рагатеег$ ипаег 
{е а{егпайпд оадб, [5 изед. Тпе аи ог5 ргемоиз депегайе4 {Теогейса!/ тодЕ! о! пе сотЬте апуе дупат/с 
5ует иййн пе поттеаг попвоопотис Йг5Е-огаег сопзвгатЕ 15 изеа 0 Ююпт ве $Еие зрасе. 

Кеуигога$: сотЬте апуе дупатис$, поп/теаг попвоопотис @г5Еогаег сопзватЕ $еерезЁ-ЦезсепЁ те од, орйта! 
аеяоп. 
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